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The  three-dimensional  microstructure  of  an  SOFC  anode  has  been  characterized  using  a  focused  ion 
beam-scanning  electron  microscope.  The  sample  preparation  and  the  experimental  milling  and  imaging 
parameters  have  been  optimized  in  order  to  obtain  a  high-quality  3D  reconstruction.  Volume  fractions 
and  the  volumetric  connectivity  of  the  individual  phases,  specific  surface  and  interface  areas  and  the 
three-phase  boundary  length  have  been  estimated.  Effective  thermal,  electronic  and  ionic  conductivities 
of  the  sample  as  well  as  the  tortuosity  of  the  solid  phases  have  been  evaluated  by  solving  the  diffusive 
transport  equation  with  an  implicit  3D  finite  difference  method. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  solid  oxide  fuel  cell  (SOFC)  technology  is  one  of  the  most 
promising  energy  conversion  systems  due  to  its  high  efficiency  and 
fuel  flexibility  [1  ].  The  materials  involved  in  the  constitution  of  such 
fuel  cells  have  become  increasingly  sophisticated,  both  in  composi¬ 
tion  and  microstructure.  One  of  the  most  common  anode  materials 
for  SOFCs  is  a  porous  Ni-YSZ  cermet.  The  anode  is  the  place  where 
the  oxidation  of  the  fuel  gas  occurs  according  to  the  reaction: 

H2  +  02“^  H20  +  2e“ 

The  most  important  requirements  for  the  anode  are: 

(1)  catalytic  activity:  the  anode  must  have  a  high  catalytic  activ¬ 
ity  for  the  oxidation  of  the  fuel  gas:  the  greater  the  number  of 
reaction  sites  the  greater  the  potential  for  chemical  reactions; 

(2)  stability:  the  anode  must  be  chemically,  morphologically,  and 
dimensionally  stable  at  the  required  operating  temperature  in 
the  fuel  gas  environment;  it  must  also  be  resistant  towards 
contaminants; 

(3)  conductivity:  the  greatest  electrical  conductivity  under  a  large 
variety  of  operating  conditions  is  desired  to  minimize  ohmic 
losses,  and  a  high  ionic  (O2-)  conductivity  is  also  expected  in 
order  to  promote  the  reaction  of  the  fuel  with  oxide  ions; 
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(4)  compatibility:  the  anode  must  be  chemically,  thermally,  and 
mechanically  compatible  with  the  other  cell  components  dur¬ 
ing  fabrication  as  well  as  under  operation; 

(5)  porosity:  the  porosity  of  the  anode  must  be  adapted  to  optimize 
mass  transport. 

For  most  of  these  conditions,  some  experimental  studies  have 
been  performed  to  estimate  the  corresponding  parameters:  vol¬ 
ume  fraction,  phase  connectivity,  specific  surface  area,  three-phase 
boundary  length  (TPBL),  tortuosity,  as  well  as  thermal,  electrical 
and  ionic  conductivities.  These  parameters  can  only  be  described 
three-dimensionally  and  thus  cannot  be  estimated  using  current 
methods  of  microstructural  analysis,  like  for  example  scanning 
electron  microscopy  (SEM)  which  only  provides  two-dimensional 
singularities  of  the  microstructure. 

Recently,  3D  reconstructions  of  SOFC  electrodes  (both  anode 
and  cathode  with  characteristic  particle  sizes  in  the  range  of 
1 00-500  nm)  have  been  performed  by  X-ray  computed  tomog¬ 
raphy  (XCT)  [2-5]  and  primarily  by  focused  ion  beam-scanning 
electron  microscopy  (FIB-SEM)  [6-22].  The  FIB  tomography 
principle  consists  in  ablating  a  structure  physically  by  FIB  slic¬ 
ing  (Ga+  ions)  followed  by  a  digital  reconstruction  based  on 
SEM  images  collected  after  each  ablation  step.  The  process 
can  be  carried  out  manually  or  automatically.  The  obtained 
images  are  then  aligned  and  superimposed  in  order  to  gen¬ 
erate  the  volume.  A  few  hundred  images  are  recorded  with 
an  inter-slice  distance  of  about  30-1 00  nm  (volume  dimension 
typically  close  to  10  x  10  x  10  p,m3).  The  voxel  dimensions  and 
the  size  of  the  analyzed  volume  must  be  optimized  according 
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Table  1 

Characteristic  data  of  the  starting  powders  and  ink. 


Components 

Suppliers 

Particle  size,  dso  (|Jim)a 

Surface  area  (m2  g-1  )a 

Viscosity  (Pas  @  10s_1)b 

TZ-8YS  powder 

Tosoh 

0.55 

6.20 

NiO  powder 

JT  Baker 

0.9 

3.30 

- 

Ink  medium 

Nextech  materials 

- 

- 

0.1 

a  Data  measured  from  Cilas  granulometer. 
b  Information  from  the  supplier. 


to  the  particle  size  and  heterogeneity  of  the  studied  sample 
[23]. 

As  described  below,  the  obtained  3D  microstructure  can  be 
either  analyzed  to  extract  some  3D  physical  parameters  or  directly 
implemented  into  numerical  models,  in  order  to  for  example  cal¬ 
culate  the  transport  properties  (of  gas,  heat,  ions  and  electrons). 
On  the  one  hand,  the  volume  data  can  be  used  to  estimate  features 
like  volume  fraction,  phase  connectivity,  interfacial  areas,  TPBL  and 
tortuosity.  The  microstructure  is  therefore  described  by  some  geo¬ 
metrical  3D  parameters.  These  parameters  are  then  included  into 
analytic  relations  originating  from  simple  models  or  experimen¬ 
tal  measurements  to  evaluate  the  electrochemical  performance  of 
the  corresponding  sample.  For  instance,  the  Tanner-Fung-Virkar 
composite  electrode  model  [24]  has  been  employed  to  predict  the 
polarization  resistance  of  SOFC  anodes  and  cathodes  [7,11].  Elec¬ 
trochemical  impedance  measurements  on  SOFC  cathodes  have  also 
been  correlated  to  3D  analysis  to  establish  an  analytical  relation 
between  the  polarization  resistance  and  the  TPBL  [14,15]. 

On  the  other  hand,  the  data  can  also  be  used  as  input  for  SOFC 
transport  models  to  assess  the  electrode  performance.  For  instance, 
the  anodic  overpotential  has  been  calculated  by  coupling  the  mass 
and  charge  transfer  and  the  electrochemical  reaction  using  the 
Lattice  Boltzmann  method  [18].  These  two  approaches  have  been 
applied  to  electrodes  with  different  YSZ/Ni  ratios  [11,12,14]  or 
subjected  to  aging  under  a  controlled  environment  including  tem¬ 
perature,  fuel  composition  or  humidity  [20,21  ]. 

The  present  study  employed  the  optimization  approach  to 
obtain  a  high-quality  3D  reconstruction  of  an  Ni-YSZ  anode 
using  FIB-SEM.  Then,  from  these  3D  data,  various  microstructural 
parameters  (volume  fraction,  phase  connectivity),  and  interfacial 
parameters  (surface  area  and  TPBL)  were  quantified.  Moreover,  by 
solving  the  diffusive  transport  equation  on  the  analyzed  volume, 
the  effective  thermal,  electrical  and  ionic  conductivities  of  the  sam¬ 
ple  and  the  tortuosity  of  the  solid  phases  could  be  estimated. 

2.  Experimental 

2A.  Preparation  of  the  anodic  cermet 

The  deposition  of  the  analyzed  NiO-8YSZ  anode  (a  three-layer 
coating  on  the  electrolyte  substrate)  was  performed  by  screen¬ 
printing  onto  circular  8  mol.%  yttria  doped-zirconia  (8YSZ)  supports 
using  ink  with  an  optimized  composition  and  viscosity.  NiO  from 
J.T.  Baker  (USA)  and  8YSZ  from  Tosoh  (Japan)  were  used  as  the  raw 
materials  to  prepare  the  NiO-8YSZ  cermet  at  varying  Ni  volume 
ratios.  Table  1  gives  the  particle  size,  surface  area  and  ink  viscosity. 

First  of  all,  in  order  to  prepare  a  composite  powder  with  an  NiO 
ratio  of  56wt.%  (corresponding  to  40vol.%  Ni),  the  starting  pow¬ 
ders  were  weighed  and  ball-milled  in  ethanol  during  18h.  Next, 
the  solvent  was  evaporated  in  an  atex  oven  during  24  h  at  90  °C. 
After  cooling,  and  to  prevent  agglomerates  from  forming,  a  200- pan 
screening  led  to  an  NiO-8YSZ  composite  powder  with  the  desired 
NiO  ratio. 

The  anode  ink  was  obtained  by  mixing  the  NiO-8YSZ  pow¬ 
der  (65.5  wt.%)  with  alpha-terpiniol  (4.9wt.%)  and  Nextech  ink 
(27.8  wt.%)  as  the  solvent  and  medium,  respectively,  and  alkyl 


ether  phosphate  acid  was  employed  as  the  dispersant  (1.8  wt.%). 
After  mixing  and  grinding  the  anode  ink  for  1 5  min,  it  was  screen- 
printed  onto  circular  zirconia  supports  through  PET1000  screens. 
The  first  layer  of  the  anode  ink  was  dried  in  an  oven  for  20  min  at 
100°C  before  a  second  layer  of  the  same  anode  composition  was 
screen-printed  onto  the  first.  The  deposition  of  the  third  layer  was 
processed  in  the  same  way.  The  electrodes  were  fired  at  1300  °C 
during  3h.  The  reduction  treatment  was  performed  at  800  °C  for 

4  h  using  a  2  °C  min-1  ramp  rate  in  Ar-Fl2  (2  vol.%,  0.5 1  min-1 ).  By 
assuming  a  complete  reduction  of  the  NiO  in  Ni,  the  8YSZ/Ni  volume 
ratio  was  expected  to  be  60/40  after  the  reduction  treatment. 

2.2.  FIB  procedure 

The  image  set  was  collected  with  a  FEI  dual-beam  Strata  400- 

5  FIB  system.  The  procedure  to  obtain  the  3D  data  from  the  initial 
sample  is  schematically  presented  in  Fig.  1  and  described  as  follows. 

2.2. 1.  Sample  preparation  before  FIB  experiment 

In  order  to  easily  distinguish  the  pores  during  SEM  observa¬ 
tion  and  to  avoid  mistakes  due  to  the  depth  of  field  in  the  SEM 
images,  the  sample  was  infiltrated  with  an  epoxy  resin  under  vac¬ 
uum  (EPOFIX,  Struers).  Subsequently,  it  was  polished  up  to  a  0.5  p,m 
finishing  and  cut  to  obtain  a  typical  thinness  of  500  p,m.  This  was 
done  to  limit  sample  drift  during  FIB  operation.  Finally  the  sam¬ 
ple  was  glued  onto  an  aluminum  stub,  and  a  conducting  gold  layer 
(~100  nm)  was  deposited  by  sputtering. 

2.2.2.  Preparation  of  the  region  of  interest  (ROI) 

The  sample  was  then  mounted  in  the  microscope  chamber  and 
adjusted  to  eucentric  height  so  that  the  exposed  surface  could  be 
imaged  directly  in  SEM  mode  without  repositioning.  To  this  end, 
as  illustrated  in  Fig.  2,  the  sample  surface  was  perpendicular  to  the 
ion  gun  and  at  an  angle  of  52°  with  respect  to  the  electron  gun. 

Prior  to  the  FIB  slicing  process,  a  representative  cross-sectional 
area  of  the  sample  was  located  by  electron  imaging  and  coated  with 
a  Pt  protective  layer  of  2-3  jxm  (Fig.  2)  using  an  in  situ  metallorganic 
ion  source.  This  was  done  in  order  to  protect  the  sample  from  acci¬ 
dental  ion  milling  and  erosion  during  ion-beam  imaging  and  also  to 
smooth  the  surface  topography.  Above  all  it  was  necessary  to  avoid, 
or  at  least  strongly  reduce  a  common  artifact  known  as  “water¬ 
falling”  or  the  “curtaining”  effect,  seen  as  an  example  on  the  cross 
section  in  Fig.  3(a)  since  it  occluded  the  entire  image.  A  Pt  deposi¬ 
tion  associated  with  sample  infiltration  and  surface  polishing  have 
thus  been  realized 

FIB  was  used  to  mill  wide  and  deep  trenches  around  the  ROI 
with  a  maximal  aperture  current  of  21  nA,  giving  rise  to  a  “U-shaped 
trench”  [  1 2  ].  A  relatively  high  dimension-free  space  around  the  ROI 
was  essential  to  prevent  the  re-deposition  of  the  sputtered  material 
during  the  sectional  milling  process.  If  the  space  around  the  ROI  was 
too  small,  images  with  high  concentrations  of  artifact  features  or 
uneven  brightness  were  obtained  (Fig.  3(b)). 

2.2.3.  Milling  and  imaging  process 

In  order  to  determine  the  depth  removed  in  the  z  direction  (see 
Fig.  2)  between  two  consecutive  sections  (interslice:  A z),  the  rel- 
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Fig.  1.  Schematic  representation  of  the  experimental  method  used  to  obtain  the  3D 
reconstruction  of  the  SOFC  anode  from  FIB  tomography. 
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Fig.  2.  SEM  image  of  the  sample  before  the  “milling  and  imaging”  procedure.  The 
stage  tilt  was  kept  the  constant  angle  of  52°  from  the  electron  beam  while  the  sample 
surface  ( xz  plane)  was  perpendicular  to  the  ion  beam  (y  direction).  The  region  of 
interest  (ROI)  is  displayed  in  the  orange  rectangle.  The  sample  was  sliced  in  the  z 
direction.  The  small  rectangle  on  the  Pt  deposition  (reference  mark)  was  used  to 
measure  the  interslice  ( Az)  and  correct  for  drifting  effects. 


The  SEM  and  FIB  operating  conditions  for  the  sequential  ion 
milling  and  electron  imaging  were  optimized.  The  FIB  milling  was 
performed  with  an  ion-beam  current  of  460  pA  at  30  kV.  For  SEM 
imaging,  a  magnification  of  12  kX  and  a  Through  the  Lens  Detector 
(TLD)  operating  in  BSE  mode  at  a  low  scan  rate  were  employed.  As 
the  difference  between  the  backscatter  coefficients  of  Ni  and  8YSZ 
were  at  a  maximum  for  an  accelerating  voltage  smaller  than  1  kV 
[25],  a  0.5-kV  accelerating  voltage  provided  an  optimal  contrast 
between  Ni  (bright),  8YSZ  (grey)  and  the  pores  (dark)  (Fig.  4(a)). 
Moreover,  this  low  voltage  rendered  possible  the  collection  of  only 
backscattered  electrons  (BSE)  emitted  very  close  to  the  imaged  sur¬ 
face. 

Indeed,  if  the  interaction  region  between  the  electron  beam  and 
the  material  is  larger  than  the  interslice  value,  the  3D  imaging  will 
not  be  correct.  Monte  Carlo  calculations  were  performed  with  the 
CASINO  software  [26]  on  each  constituent  of  the  studied  electrode: 
Ni,  8YSZ  and  resin  (density:  8.9,  5.9  and  1.1  gem-3  respectively) 
and  displayed  that  for  a  0.5-kV  accelerating  voltage,  BSEs  arose 
from  a  maximum  depth  of  6  nm,  which  was  much  smaller  than 
the  measured  interslice  value  equal  to  lOOnm.  In  this  case,  115 
images  were  recorded  using  a  manual  procedure  along  the  sample 
thickness  (Fig.  4(a)).  Including  the  in  situ  sample  preparation  and 
sectioning,  the  total  processing  time  was  approximately  18  h. 


ative  recession  of  the  depth-profiled  surface  was  measured  and 
compared  to  the  original  cross-section  based  on  a  reference  mark 
(rectangle)  on  the  upper  part  of  the  Pt  layer.  Since  the  “milling  and 
imaging”  procedure  was  very  time-consuming,  uncorrected  drift  of 
electron  beam,  stage  and  sample  could  greatly  affect  the  interslice 
real  value.  For  this  reason  a  correction  procedure  was  utilized  [23]. 
Drift  values  in  the  x  and  z  directions  were  measured  and  corrected 
after  each  slicing  and  imaging  loop.  For  this  purpose,  an  image  con¬ 
taining  the  reference  mark  (rectangle  in  Fig.  2)  was  acquired  using 
the  ion-induced  secondary  electron  detector  (CDEM:  continuous 
dynode  electron  multiplier)  located  perpendicularly  to  the  sam¬ 
ple  surface.  Before  recording  an  image  used  for  the  reconstruction, 
a  CDEM  image  was  registered  and  the  ion  beam  was  eventually 
deflected  to  allow  perfect  matching  between  the  actual  and  original 
positions  of  the  milled  rectangle. 


2.2.4.  Data  processing 

Data  processing  includes  the  following  steps:  (i)  alignment  of 
the  consecutives  slices,  (ii)  correction  of  the  dimensions  taking  into 
account  the  tilt  of  the  electron  beam,  (iii)  thresholding  of  the  grey 
levels  and  labeling  of  phases,  (iv)  resampling  of  the  data  to  obtain 
cubic  voxels  and  (v)  3D  image  generation.  Most  of  the  processing 
steps  were  performed  using  the  AVIZO  6.2  software  [27].  Details  of 
each  step  are  given  below. 

The  still  remaining  drift  components  of  the  electron  beam  in  the 
x  and  y  directions  were  corrected  by  applying  least  square  fitting 
algorithms  to  achieve  a  correct  alignment  of  the  images.  Finally, 
the  resulting  volume  was  cropped  to  keep  only  the  ROI  common  to 
all  slices.  The  absolute  dimensions  in  the  x  direction  were  obtained 
from  a  calibrated  SEM  magnification.  Because  of  oblique  SEM  imag¬ 
ing  at  an  angle  of  52°,  distances  in  they  direction  were  corrected  for 
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Fig.  3.  SEM  images  of  a  Ni-8YSZ  anode  illustrating  some  common  FIB  artifacts:  (a)  “curtaining”  effect  due  to  a  too  small  (<1  (Jim)  Pt  thickness,  (b)  deposition  of  the  sputtered 
material  (arrows)  and  shadowing  effect  (uneven  brightness  from  the  top  to  the  bottom  of  the  image)  occurring  when  trenches  around  the  ROI  are  not  sufficiently  deep  and 
wide. 


projection  effects  using yreal  =ysEM_image/sin  (52°).  In  certain  cases, 
a  slight  tilting  occurred  because  of  imperfect  straight  milling  in  the 
XZ  plane  (see  Fig.  2);  however,  since  the  460  pA  milling  current  was 
small,  this  effect  was  neglected.  The  3D  image  dimensions  were 
finally  8.66  x  9.79  x  1 1 .41  pan3  (volume  =  967.36  pan3 ). 

Even  though  a  good  and  reproducible  contrast  between  the 
three  phases  could  be  obtained,  the  boundaries  between  them  were 
generally  not  perfectly  sharp  lines,  but  rather  gradual  transition 
zones  a  few  voxels  wide.  Consequently,  a  median  filter  and  a  delin¬ 
eating  operation  were  applied  to  the  entire  stack  of  images  to  obtain 
a  good  separation  between  each  phase.  Segmentation  of  the  data 
(assigning  a  phase  to  a  voxel)  was  then  realized  by  grey  level-based 
thresholding  (Fig.  4(b)). 

The  initial  dimension  of  voxels  was  10.4  x  13.2  x  100  nm3.  How¬ 
ever,  the  dataset  was  resampled  to  obtain  approximately  one  billion 
10  x  10  x  10nm3 cubic  voxels  using  a  cubic  interpolation  method. 
This  was  due  to  the  fact  that  the  following  analysis  method  (TPBL 
and  tortuosity  measurements)  required  cubic  voxels.  The  resam¬ 
pling  step  was  carried  out  after  the  segmentation  of  the  data  to 
avoid  artifacts  caused  if  this  step  was  applied  directly  to  the  grey 
level  data. 

3.  Quantification  results  and  discussion 

First,  volume  fractions  of  each  phase  and  their  corresponding 
connectivity  across  the  analyzed  volume  were  determined.  The  sur¬ 


face  area  of  interfacial  planes  between  two  phases  as  well  as  the 
lengths  of  the  three-phase  boundary  lines  were  evaluated.  Then, 
the  effective  thermal,  electronic  and  ionic  conductivities  of  the 
electrode  and  the  tortuosity  of  the  solid  phases  were  estimated  by 
numerically  solving  the  diffusive  transport  equation.  Most  of  this 
quantitative  analysis  was  performed  with  laboratory-developed  C 
programs  on  a  personal  computer  [dual  CPU  Intel®  Xeon®  E5440  @ 
2.83  GHz,  32  GB  RAM]. 

3.1.  Volume  fraction 

The  volume  fractions  calculated  from  3D  data  provided  the  fol¬ 
lowing  results:  41  %  for  pores,  33%  for  8YSZ  and  26%  for  Ni.  The  solid 
volume  fractions  have  been  estimated  to  44%  Ni  and  56%  8YSZ  in 
the  studied  anode,  which  agreed  well  with  the  expected  values  (i.e., 
40  vol.%  Ni/60  vol.%  8YSZ).  This  result  provided  evidence  of  the  good 
representativity  of  the  analyzed  volume. 

3.2.  Volumetric  connectivity 

In  order  to  verify  the  connectivity  of  each  phase,  a  laboratory- 
developed  program  based  on  a  Hoshen-Kopelman  algorithm  was 
employed  [28].  Voxel  after  voxel,  the  program  verified  the  connec¬ 
tivity.  The  voxels  that  were  linked  together  by  a  face  were  judged  to 
be  connected  to  form  a  cluster.  When  two  voxels  were  only  linked 
by  a  vertex  or  an  edge,  they  were  considered  not  to  be  connected. 


Fig.  4.  Images  of:  (a)  one  of  the  original  115  serial-sections  collected  with  a  BSE  detector,  so  that  Ni  appears  brighter  than  8YSZ  and  the  impregnated  pores  appear  dark,  (b) 
final  3D  image  after  alignment,  delineating,  cropping  and  adjustment  of  voxel  size  operations  showing  Ni  (bright),  8YSZ  (grey)  and  pore  (dark)  phases. 
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Fig.  5.  An  Ni  cluster  representation  displaying  that  the  main  part  (87.4%)  of  the  Ni 
phase  was  percolated,  i.e.,  in  contact  with  the  six  boundary  faces  of  the  volume  (in 
yellow). 


This  cluster  neighborhood  rule  is  commonly  used  in  connectivity 
analysis  [29].  Finally,  if  a  cluster  was  connected  to  the  six  boundary 
faces  of  the  studied  volume,  it  was  defined  as  “percolated”,  while 
others  remained  isolated  (Fig.  5). 

The  performed  analysis  revealed  that  the  three  phases  were 
mainly  constituted  of  a  very  large  “percolated”  cluster  which  rep¬ 
resented  99.8%,  99.1%  and  87.4%  of  the  pore,  8YSZ  and  Ni  volumes, 
respectively.  This  significant  percolation  of  all  three  phases  sug¬ 
gested  a  good  ability  for  the  transport  of  gases  in  the  pores,  as  well 
as  of  charge  carriers  in  the  solid  phases. 


3.3.  Interface  properties 

Electrochemical  and  catalytic  activities  of  an  electrode  are  gen¬ 
erally  described  by  a  single  parameter:  the  TPBL  which  is  correlated 
to  the  number  of  regions  where  the  electrochemical  reactions  can 
take  place.  Flowever,  the  oxidation  of  the  fuel  involves  numerous 
reaction  steps  which  are  not  necessarily  located  at  the  TPB  [20]. 
For  example,  adsorption  and  dissociation  of  hydrogen  occur  on  the 
nickel  surface  because  of  its  good  catalytic  properties.  Thereby,  only 
the  Ni-pore  interface,  as  opposed  to  the  entire  nickel  surface,  is 
exposed  to  the  gas  and  can  therefore  react  with  the  fuel.  Flowever, 
not  all  of  the  interfaces  may  be  active  due  to  non-contiguous  regions 
(isolated  cluster).  Consequently,  only  the  specific  surface  and  inter¬ 
face  areas  inside  the  percolated  volumes  have  been  calculated,  and 
the  TPBl  was  estimated  for  both  the  total  and  percolated  volumes 
for  comparison. 

The  specific  surface  areas  (surface  area  normalized  to  the  sam¬ 
ple  volume)  were  4.27  |xm2  |xm-3  and  4.24  p,m2  pan-3  for  the 
porous  and  8YSZ  phases,  respectively,  and  2.33  pm2  p,nrr3  for 
the  Ni  phase.  The  pie  diagram  of  Fig.  6  indicates  the  specific 
interface  areas  (interface  area  normalized  to  the  sample  vol¬ 
ume).  For  instance,  the  specific  surface  area  of  the  Ni  phase 
(2.33  pm2  p,nrr3)  was  divided  into  Ni-pore  (1.18  pm2  pm-3)  and 
Ni-8YSZ  (1.15  pm2  pm-3)  interfaces.  Thus,  about  50%  of  the  total 
nickel  surface  was  exposed  to  the  porous  phase  and  could  be  used 
for  surface  catalytic  reactions  with  the  fuel. 


Interface 

NI-8YS2 


Interface 

Ni-Pore 


Fig.  6.  Pie  diagram  representing  the  specific  surface  areas  of  pores,  8YSZ  and  Ni, 
and  the  specific  interface  areas  between  the  three  neighboring  phases  (surface  and 
interface  are  normalized  by  volume  sample).  The  active  TPBL  is  7.4  |xm  ijutt3. 


The  TPBs  are  the  regions  where  the  three  phases  (pores,  8YSZ, 
Ni)  met,  thus  corresponding  -  on  the  3D  image  -  to  the  edges  of 
voxels  common  to  the  three  phases  [4].  A  program  has  been  devel¬ 
oped  to  identify  and  count  these  edges.  The  following  procedure 
was  applied  to  each  voxel  by  scanning  the  volume  in  the  positive 
direction  (+x,  +y,  +z),  from  the  “first”  voxel  (1,1,1)  to  the  last  one 

(^max>  Ymax*  Zmax)- 

Considering  the  schematic  representation  of  Fig.  7(a),  if  in  the 
xy  plane  at  least  two  of  the  three  neighboring  voxels  have  a  nature 
differing  from  that  of  the  studied  voxel,  the  corresponding  edge  of 
contact  between  the  three  phases  (designated  edge  1  in  Fig.  7(a)) 
is  labeled  a  TPB.  The  same  method  was  applied  in  the  xz  plane  and 
yz  planes,  leading  respectively  to  edge  2  and  edge  3  (Fig.  7(a)).  The 
TPBl  was  calculated  by  summing  edge  1,  edge  2  and  edge  3  all 
over  the  dataset,  and  this  value  was  then  normalized  to  the  studied 
volume.  By  using  this  procedure,  double  counting  of  a  TPB  interface 
unit  could  be  avoided  [4]. 

This  TPB  analysis  procedure  was  first  applied  to  the  initial 
volume  and  yielded  a  TPB  density  of  1 1 .2  p,m  [xrrr3.  The  3D  visual¬ 
ization  of  these  TPBs  was  not  easy  since  elementary  voxels  have  to 
be  represented  into  edges.  Fig.  7(b)  represents  the  voxels  assigned 
to  pores  containing  at  least  one  TPB,  because  the  graphic  software 
does  not  support  edge  voxel  representation.  Flowever,  more  than 


Fig.  7.  (a)  Determination  of  TPBL  from  the  edge  segment  length  measurement:  in 
the  xy  plane,  if  at  least  two  of  the  three  neighboring  voxels  correspond  to  a  differ¬ 
ent  phase  than  the  considered  voxel  (in  yellow),  the  corresponding  edge  of  contact 
between  the  three  phases  (edge  1)  is  labeled  as  a  TPB.  The  same  method  is  applied 
in  the  xz  plane  and  yz  planes,  giving  edges  2  and  3,  respectively.  Finally,  edges  1, 
2  and  3  are  summed  and  divided  by  the  volume  dimension  to  calculate  TPBL.  (b) 
A  3D  map  of  the  “active”  TPBs:  all  voxels  assigned  to  pores  and  containing  at  least 
one  active  TPB  are  represented.  (For  interpretation  of  the  references  to  color  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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one  TPB  can  belong  to  a  pore  voxel  and  as  a  consequence  it  induces 
an  under-representation  of  the  total  number  of  TPBs  (statistically 
for  100  TPB  edges,  only  60  voxels  are  represented). 

The  “activity”  of  a  TPB  is  subject  to  numerous  definitions 
in  the  literature.  A  majority  of  authors  have  defined  a  TPB  as 
“active”  if  it  comprises  three  percolated  phases  (pore,  Ni,  8YSZ) 
[4,6,8,10,11,16,17,20].  A  TPB  can  also  be  defined  as  “active”  if  the 
three  following  conditions  are  verified:  (1)  there  is  a  contiguous 
network  of  8YSZ  from  the  TPB  of  interest  to  the  electrolyte,  (2 )  there 
is  a  contiguous  network  of  Ni  from  the  TPB  of  interest  to  the  current 
collector  and  (3)  there  is  a  contiguous  network  of  pores  from  the 
TPB  of  interest  to  the  inlet  of  the  gas  [30].  However,  this  last  defini¬ 
tion  can  only  be  used  if  the  reconstructed  volume  is  extended  from 
the  electrolyte  to  the  collector. 

The  first  definition  was  here  applied  to  the  TPB  research  proce¬ 
dure  on  the  percolated  volumes  as  mentioned  before.  The  “active” 
TPBl  was  calculated  and  found  equal  to  7.4  p,m  |jinrr3,  which  rep¬ 
resents  66%  of  the  TPBl  calculated  before  and  corresponding  to 
the  total  volume.  This  value  was  relatively  high  compared  to  oth¬ 
ers  reported  in  the  literature  (based  on  the  same  TPBL  calculation 
method)  which  were  in  the  range:  2-1 0  p,m  p,m-3  [6-8,1 5,1 6].  This 
suggests  an  elevated  electrochemical  performance  for  the  studied 
electrode.  The  calculation  of  the  3D  parameters:  volume  fraction, 
connectivity,  and  interface  parameters  were  done  simultaneously 
and  completed  in  approximately  5  min. 

It  is  worth  noticing  that  certain  mathematical  corrections  [11] 
or  other  TPB  research  approaches  like  the  “centroid  method” 
or  the  “voxel  expansion  method”  [17,18]  have  been  reported  in 
the  literature  in  order  to  take  into  account  the  “step-like”  pat¬ 
tern  of  voxels  edges,  giving  more  accurate  results.  The  main 
point  is  to  compare  tortuosity  values  evaluated  by  the  same 
method. 


3.4.  Effective  conductivities  and  tortuosity  calculations 


Diffusion  through  a  medium  is  usually  described  by  Eq.  (1 ): 

J  =  -D0VO  (1) 

where  V0  is  the  gradient  of  (i)  the  temperature  in  the  case  of  heat 
conduction,  (ii)  the  electronic  potential  in  the  case  of  electronic  con¬ 
duction,  (iii)  the  ionic  potential  in  the  case  of  ionic  conduction  and 
(iv)  the  concentration  in  the  case  of  mass  diffusion,  j  is  the  flux  of 
heat,  electron,  ions  or  gas  and  D0  is  (i)  the  thermal  conductivity  for 
heat  conduction,  (ii)  the  electronic  conductivity  for  electron  con¬ 
duction,  (iii)  the  ionic  conductivity  for  ionic  conduction  transport 
and  (iv)  the  diffusion  coefficient  for  the  mass  diffusion. 

In  the  case  of  a  multiphase  material,  Eq.  (1 )  was  empirically  cor¬ 
rected  [31,34,40]  by  the  diffusivity  factor  (&)  to  take  into  account 
the  presence  of  the  phases  in  which  the  transport  does  not  occur 
(for  example,  the  solid  phase  in  the  case  of  the  mass  transport).  Eq. 
(1)  becomes: 

j  =  -D  V0  (2) 

where  D  is  defined  as  the  diffusion  coefficient  in  the  material  with 
D  =  'PD0.  The  diffusivity  factor  (^)  is  linked  to  the  tortuosity  (r)  as 
follow  [40]: 


here,  e  represents  the  phase  fraction  in  which  the  diffusion  takes 
place. 

Since  D0  corresponds  to  the  thermal  conductivity,  as  well  as 
the  electronic  conductivity,  the  ionic  conductivity  and  the  diffu¬ 
sion  coefficient,  its  values  can  be  obtained  from  the  literature.  To 
calculate  D,  as  mentioned  above,  the  homogenization  approach  was 
performed  through  the  determination  of  an  effective  parameter 
Deff,  similar  to  D  but  for  the  equivalent  homogeneous  medium.  In 
this  approach,  the  ballistic  behavior  was  neglected  since  the  mean 
free  path  of  the  diffusing  species  was  assumed  smaller  than  the 
characteristic  size  of  the  diffusing  medium. 


3.4A.  Definition  of  tortuosity 

The  term  tortuosity  is  typically  introduced  as  some  kind  of 
“fudge”  factor  in  macroscopic  transport  equations  to  take  into 
account  the  complex  transport  paths  in  a  heterogeneous  media. 
In  a  geometrical  approach,  the  tortuosity  (r)  is  defined  as  the  ratio 
of  the  “effective  average  path”  (Leff)  of  a  fluid  particle  to  the  cor¬ 
responding  straight  and  shortest  distance  L  along  the  direction  of 
macroscopic  flux  [31  ].  Tortuosity  is  a  parameter  generally  used  to 
describe  mass  transport  phenomena  inside  a  porous  phase,  but  it 
can  also  be  extended  to  other  transport  processes  including  elec¬ 
tronic,  and  ionic  in  the  solid  phases.  Some  attempts  to  measure 
this  3D  parameter  have  already  been  made  experimentally  through 
a  gas  counter  diffusion  method  [32]  or  saturation  current  den¬ 
sity  measurement  [33].  However,  the  most  reported  procedure  for 
estimating  the  tortuosity  is  based  on  the  3D  analysis  of  a  porous 
medium  from  X-ray  or  FIB  tomography  [4,6,34]. 

In  most  cases,  tortuosity  values  can  be  deduced  from  a 
homogenization  approach.  In  the  case  of  a  heterogeneous  (real) 
material  with  anisotropic  properties  and  geometrical  complexity, 
the  homogenization  technique  is  an  efficient  method  to  obtain 
averaged  properties  equivalent  to  those  of  a  homogeneous  medium 
[35,36].  This  methodology  has  been  reported  for  numerical  recon¬ 
structed  composite  electrodes  (randomly  distributed  spherical 
particles)  [6,37,38].  An  alternative  approach  is  to  apply  a  Lattice 
Boltzmann  Method  (LBM)  to  solve  a  diffusion  equation  and  calcu¬ 
late  a  phase  tortuosity.  Iwai  et  al.  [17]  compared  a  random-walk 
method  with  an  LBM  approach  and  similar  tortuosity  values  were 
obtained.  In  the  present  work,  the  homogenization  approach  was 
used  by  solving  the  diffusive  transport  equation  by  means  of  a  finite 
difference  method  [39]. 


3.4.2.  Effective  diffusion  coefficient  calculation 

In  the  absence  of  external  sources  and  in  the  case  of  steady- 
state  conditions,  the  distribution  of  0  (for  all  conduction  or  diffusion 
phenomena)  in  the  material  is  governed  by  the  diffusive  transport 
equation  [39]: 

V-(-DV0)  =  0  (4) 

This  equation  was  first  numerically  solved  using  an  implicit 
finite  difference  method  [6,39]  to  obtain  the  distribution  of  tem¬ 
perature,  potential  (electronic  and  ionic)  and  concentration.  The 
following  boundary  conditions  were  fixed:  for  x*  =  0,  0  =  0\  and  for 
Xi=*i max  0  =  02\  where  x*  corresponds  to  either  x,  y  or  z  (along  the 
volume  axes,  cf.  Figs.  2  and  4(b)  for  the  direction  of  axes)  while 
V0  =  6  on  the  other  faces. 

Once  the  0  distribution  within  the  volume  was  obtained,  the 
average  conductive  flux  (0)  over  a  cross-section  of  the  surface  area 
S  normal  to  the  flux  (in  the  direction  of  the  unit  vector  n)  could  be 
calculated  according  to: 

<t>  =  (f  -DmxelV6  hdS  (5) 

where  Dvoxei  is  the  diffusion  coefficient  characteristic  of  the  phase 
(depending  of  the  considered  phenomenon:  heat,  electronic  or 
ionic  conduction  and  mass  diffusion)  present  in  the  considered 
voxel  (Table  2). 

For  the  equivalent  homogeneous  material,  the  average  conduc¬ 
tion  flux  0  can  be  written  [39]: 

0  =  DeffS^-Al  (6) 

max 
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Table  2 

Thermal,  ionic  and  electronic  conductivities  at  1123  K  of  the  pores,  8YSZ  and  Ni 
[48-52]. 


Pores 


8YSZ 


Ni 


MWm-1  K"1) 
o  ionic  (10-3  Scm-1) 
a  electronic  (S  cm-1 ) 


(H2)  0.48  [48] 

0 

0 


2  [49] 

14  [51  ] 
10-6  [52] 


50  [50] 
10-5  [52] 
104  [52] 


Before  calculating  Deff  in  the  x*  direction,  it  is  important  to 
ensure  that  the  solution  has  converged  properly.  In  this  aim,  the 
value  of  0  was  evaluated  on  each  slice  (in  the  x,  direction)  and  a 
mismatch  of  10-3  between  the  smallest  and  the  largest  0  value 
was  chosen  as  the  convergence  criteria. 

Since  the  boundary  conditions  were  identical  for  the  equivalent 
homogeneous  material  and  the  real  material,  the  fluxes  0  must  be 
equal  and  combining  Eqs.  (5)  and  (6)  gives  Deff  in  the  x*  direction: 


Deff  = 


S(0,-02) 


IL 


^voxel  ^  1  dS 


(7) 


In  the  studied  anode,  the  typical  characteristic  size  of  pore,  Ni 
and  8YSZ  particles  was  in  the  micrometer  range  as  approximated 
from  the  3D  images  of  Fig.  4(b).  In  the  case  of  electron,  ion  and 
heat  transport,  the  corresponding  mean  free  path  was  smaller  than 
lOOnm  [41-43]  and  the  bulk  diffusion  regime  assumption  could 
be  validated.  However,  this  assumption  was  no  longer  available  for 
H2  diffusion  since  the  mean  free  path  of  H2  molecules  in  a  porous 
medium  (Zh2  )  is  too  large. 


Zh2  can  be  approximated  by  the  following  expression  [40]: 


/h2  = 


RT 


(8) 


where  Na  is  Avogadro’s  number,  R  is  the  gas  constant,  T  is  the 
temperature,  dH2  is  the  diameter  of  the  H2  molecule  and  P  is  the 
pressure. 

At  an  operating  temperature  of  1123K  and  a  pressure  of  105 
Pa,  corresponding  to  the  electrochemical  test  and  by  assuming  that 
dn2  =  0-2  nm,  the  mean  free  path  of  hydrogen  molecules  was  close 
to  Zh2  =  900  nm.  This  value  is  comparable  to  the  characteristic  size 
of  the  pores.  Consequently,  H2  diffusion  in  the  studied  anode  was 
affected  by  the  Knudsen  effect  which  takes  into  account  collisions 
of  gaseous  molecules  with  the  solid  phases.  In  order  to  correctly 
estimate  the  tortuosity  of  the  porous  phase,  it  was  necessary  to 
use  a  model  which  included  both  bulk  and  Knudsen  diffusion  pro¬ 
cesses,  like  the  Dusty  gas  model  [44].  Consequently,  only  heat,  ion 
and  electron  diffusive  transports  were  calculated  in  this  work.  The 
results  are  given  in  Table  3  where  the  different  values  of  Deff  are 
reported  under  the  form  of  thermal,  electronic  and  ionic  effective 
conductivities.  The  entire  calculation  was  completed  in  approxi¬ 
mately  7  h. 

The  effective  thermal  conductivity  appeared  quite  isotropic 
with  a  slightly  lower  value  in  the  z  direction,  close  to  4  Wm-1 I<-1 . 
This  value  was  within  the  same  order  of  magnitude  as  those  gener¬ 
ally  reported  for  SOFC  anodes  [45-47].  As  observed  in  Table  3,  the 
effective  electronic  conductivity  appeared  to  be  anisotropic:  the 
(^electronic  value  in  the  z  direction  was  smaller  than  the  ones  in  the 
x  and  y  directions.  This  could  be  explained  by  either  a  dissimilar  Ni 


Table  3 

Thermal,  electronic,  and  ionic  effective  conductivities  in  the  x,  y  and  z  directions 
from  the  calculation  of  Deff  at  T=  1 123 1<. 


X 

y 

z 

ATWm-MC1) 

4.23 

4.54 

3.27 

o  ionic  (10-3  Scm-1) 

1.43 

1.19 

1.03 

o  electronic  (S  cm-1 ) 

281 

252 

66 

x  distance  (pm) 


y  distance  (|jm) 


Fig.  8.  The  Ni  volume  ratio  in  the  x,  y  and  z  directions  of  the  studied  volume. 


concentration  along  the  z  direction  as  compared  to  along  the  x  and 
y  directions,  or  by  a  different  distribution  of  Ni  in  the  successive 
planes  perpendicular  to  the  z  axis. 

In  the  graph  of  Fig.  8,  the  Ni  ratio  in  all  the  resampling  planes 
is  plotted  in  the  three  directions,  showing  few  variations  between 
them.  As  a  consequence,  the  Ni  distribution  plays  the  most  impor¬ 
tant  role  for  the  anisotropic  electrical  behavior.  Experimental 
four-point  electrical  measurements  performed  at  1123K  on  the 
anode  surface  led  to  a  value  equal  to  81  Scm-1.  This  value  was 
within  the  same  order  of  magnitude  as  its  calculated  counterparts 
and  the  difference  between  them  can  have  many  explanations,  both 
from  an  experimental  point  of  view  and  with  regard  to  the  signi¬ 
fication  of  an  average  value  compared  to  those  given  in  the  three 
directions. 

The  ionic  conductivity  in  the  studied  anode  was  quite  isotropic 
and  close  to  10-3Scm-1.  Since  a  value  of  1.4  x  10-2Scm-1  had 
been  used  for  the  calculation,  the  microstructural  effect  for  the  loss 
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Table  4 

Tortuosity  of  the  8YSZ  and  Ni  phases  in  the  x,  y  and  z  directions  calculated  by  solving 
the  transport  equation  by  the  finite  difference  method. 


X 

y 

z 

T8YSZ 

1.79 

1.96 

2.10 

rNi 

3.04 

3.18 

6.24 

of  conductivity  was  near  a  factor  10.  This  factor  might  be  different 
for  another  8YSZ  content  and  distribution  inside  the  material. 


of  tortuosity  able  to  provide  an  interesting  comparison  between 
materials  with  different  Ni  ratios,  different  sintering  temperature, 
fabricated  by  different  routes  or  after  aging. 

In  a  forthcoming  paper,  FIB  tomography  experiment  and  cal¬ 
culation  procedures  will  be  presented  for  Ni-8YSZ  cermets  of 
various  compositions,  ranging  from  70/30  to  55/45  vol.%  of  8YSZ/Ni. 
Microstructural  parameters,  interfacial  features  and  tortuosity  of 
the  solid  phases  will  be  discussed  as  a  function  of  the  composition 
of  the  electrode. 


3.4.3.  Tortuosity  results 

As  seen  in  Table  2,  the  ionic  conductivity  of  Ni  and  pores,  and 
the  electronic  conductivity  of  8YSZ  and  pores  were  negligible.  It 
was  thus  possible  to  deduce  the  tortuosity  of  the  solid  phases  from 
the  effective  electrical  and  ionic  conductivity  values  by  combining 
Eqs.  (1  )-(3).  Eqs.  (9)  and  (10)  allow  the  calculation  of  the  tortuosity 
of  the  8YSZ  and  Ni  phases  according  to: 


r8YSZ  = 


flo(8YSZ) 

Deff 


rNi  = 


^O(Ni) 

Deff 


(9) 

(10) 


where  £8ysz  and  £Ni  are  the  8YSZ  and  Ni  volume  fractions,  respec¬ 
tively  (cf.  Section  3.1),  D0  (8YSZ)  and  D0  (Ni)  are  the  electrical  and 
ionic  conductivities  of  the  8YSZ  and  Ni  phases  (cf.  Table  2).  More¬ 
over,  r8Ysz  and  rNi  can  be  calculated  in  the  x,y  and  z  directions  (see 
Table  4). 

The  tortuosity  of  the  Ni  phase,  considered  as  the  only  electronic 
conductor,  was  higher  in  the  z  direction  (6.24)  as  compared  to  in 
the  x  and  y  ones  (3.04  and  3.18,  respectively).  This  was  not  surpris¬ 
ing  since  the  value  was  directly  deduced  from  Deff.  Using  the  same 
considerations,  the  tortuosity  of  the  8YSZ  phase,  considered  as  the 
only  ionic  conductor,  gave  more  isotropic  values  as  compared  to 
other  reported  ones  [12]. 


4.  Conclusions 


This  work  has  focused  on  the  3D  analysis  of  the  microstruc¬ 
ture  in  an  Ni-8YSZ  SOFC  anode  reconstructed  by  FIB  tomography. 
Sample  preparation,  milling  and  imaging  conditions  as  well  as 
data  processing  have  been  optimized  to  obtain  a  high-quality 
representative  3D  microstructure  with  the  following  dimen¬ 
sions:  8.66  x  9.79  x  11.41  |jLm3.  Laboratory-developed  programs 
rendered  possible  the  calculation  of  some  microstructural  parame¬ 
ters  such  as  the  volume  fraction,  connectivity,  specific  surface  and 
interface  areas  and  TPB  length  including  the  proportion  of  active 
and  non-active  ones.  These  calculations  demonstrated  that  for  an 
anode  with  an  initial  8YSZ/NiO  volume  ratio  of  60/40,  which  after 
reduction  gave  volume  relative  amounts  of  41%  of  pores,  33%  of 
8YSZ  and  26%  of  Ni,  all  the  phases  were  highly  percolated.  The  active 
TPB  length  was  estimated  to  7.4  [xm  [xm-3  which,  in  comparison  to 
results  reported  in  the  literature,  corresponds  to  an  effective  anode. 

By  solving  the  diffusive  transport  equation  with  finite  difference 
calculations,  the  effective  thermal,  electronic  and  ionic  conductiv¬ 
ities  could  be  determined.  The  tortuosity  of  the  solid  phases  was 
deduced  from  these  calculations.  Flowever,  the  tortuosity  of  the 
porous  phase  cannot  be  easily  obtained  from  a  “pure”  diffusive 
mechanism  due  to  the  length  of  the  mean  free  path  of  the  fuel  gas, 
which  leads  to  Knudsen  diffusion,  is  difficult  to  take  into  account. 
For  this  reason,  as  in  the  finite  difference  calculation  presented 
in  this  paper,  the  tortuosity  depends  on  the  nature  of  the  species 
that  actually  diffuse  in  the  considered  medium.  No  results  have 
been  given  for  the  tortuosity  of  the  porous  phase.  Some  models  are 
evaluated  to  eventually  propose  a  simple  geometrical  description 
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